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We analyze gap solitons in trapped Bose-Einstein condensates (BECs) in optical lattice potentials 
under Feshbach resonance management. Starting with an averaged Gross-Pitaevsky (GP) equation 
with a periodic potential, we employ an envelope wave approximation to derive coupled-mode equa- 
tions describing the slow BEG dynamics in the first spectral gap of the optical lattice. We construct 
exact analytical formulas describing gap soliton solutions and examine their spectral stability using 
the Ghebyshev interpolation method. We show that these gap solitons are unstable far from the 
threshold of local bifurcation and that the instability results in the distortion of their shape. We 
also predict the threshold of the power of gap solitons near the local bifurcation limit. 

PACS numbers: 03.75.Lm, 03.75.Nt, 05.45.-a 



INTRODUCTION 

At sufficiently low temperatures, particles in a dilute 
boson gas can condense in the ground state, forming a 
Bose-Einstein condensate (BEC) 0. Under the typical 
confining conditions of experimental settings, BECs are 
inhomogeneous and the number of condensed atoms (iV) 
ranges from several thousand (or less) to several million 
(or more). The magnetic traps that confine them are 
usually approximated by harmonic potentials. There are 
two characteristic length scales: the harmonic oscillator 
length ttho = \/f^l ("niLOho) [which is on the order of a 
few microns], where LUho = {i^xi^yi^zY^^ is the geometric 
mean of the trapping frequencies, and the mean healing 
length X = 1/ \/87r|a|n [which is on the order of a mi- 
cron], where n is the mean particle density and a, the 
(two-body) s-wave scattering length, is determined by 
the atomic species of the condensate. Interactions be- 
tween atoms are repulsive when a > and attractive 
when a < 0. For a dilute ideal gas, a « 0. 

If considering only two-body, mean-field interactions, 
a dilute Bose-Einstein gas can be modeled using a cu- 
bic nonlinear Schrodinger equation (NLS) with an exter- 
nal potential; this is also known as the Gross-Pitaevsky 
(GP) equation. BECs are modeled in the quasi-one- 
dimensional (quasi-lD) regime when the transverse di- 
mensions of the condensate are on the order of its heal- 
ing length and its longitudinal dimension is much larger 
than its transverse ones The GP equation for the 
condensate wavefunction ip{x,t) takes the form 
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(1) 



where \ip\'^ is the number density, V{x) is the external 
trapping potential, g = [47r/i^a/TO] [1 -I- C'(C^)] is pro- 
portional to the two-body scattering length, and C = 



VlV'l^lal^ is the dilute gas parameter 

m 

Experimentally realizable potentials V{x) include har- 
monic traps, quartic double-well traps, optical lattices 
and superlattices, and superpositions of lattices or super- 
lattices with harmonic traps. The existence of quasi- ID 
( "cigar-shaped" ) BECs motivates the study of lower di- 
mensional models such as Eq. (^. We focus here on the 
case of spatially periodic potentials without a confining 
trap along the dimension of the lattice, as that is of par- 
ticular theoretical and experimental interest. For exam- 
ple, such potentials have been used to study Josephson 
effects 0, squeezed states 0, Landau-Zener tunneling 
and Bloch oscillations [3, period- multiplied wavefunc- 
tions Q , and the transition between superfluidity and 
Mott insulation at both the classical ^.8] and quantum 
levels. Moreover, with each lattice site occupied by one 
alkali atom in its ground state, a BEC in an optical lat- 
tice shows promise as a register in a quantum computer 

113. 

The properties of BECs — including their shape, col- 
lective excitations, statistical fluctuations, and the for- 
mation and dynamics of their solitons and vortices — are 
determined by the strength and sign of their two-body 
atomic interactions a. This scattering length, and hence 
the coefficient of the nonlinearity in the GP equation, 
can be adjusted in both sign and magnitude (over a 
large range) by minutely adjusting a magnetic field in 
the vicinity of a so-called "Feshbach resonance" 0, ^3 ■ 

A Feshbach resonance is an enhancement in the scat- 
tering amplitude of a particle incident on a target when 
the energy of the former is approximately that needed to 
create a quasibound state of the two-particle system. If 
a pair of ultracold atoms has a molecular bound state 
near zero energy, then during collisions they stick to- 
gether for a little while as they undergo a Feshbach res- 
onance. While few molecules have bound states at such 
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energies, one can adjust the relative energies of atoms and 
molecules with different magnetic moments by applying 
a magnetic field. With such "Zeeman tuning", one can 
move the atomic energy from just above the resonance to 
just below it, so that the scattering length diverges and 
changes sign from positive to negative across the reso- 
nance. 

As a result of the control this procedure gives over con- 
densate properties, the manipulation of ultracold atoms 
using Feshbach resonances has become among the most 
active research areas in experimental atomic physics. 
Feshbach resonances have provided a key for creating 
molecular BECs, generating solitons and atom-molecule 
coherence, stabilizing or destabilizing BECs, and creating 
novel Fermi liquids [isl 113. Iislll^. For example, it was 
recently shown that near a Feshbach resonance, a quan- 
tum phase transition occurs between a regime with both 
atomic and molecular condensates and one with only 
molecular condensates 0. As pointed out in Ref. [l8| . 
this transition should be much easier to observe for con- 
densates loaded into optical lattice potentials. 

In Feshbach resonance management, which was moti- 
vated by similar techniques in fiber optics the BEC 
scattering length is varied periodically in time. This 
yields dynamically interesting soliton solutions, such as 
breathers jl^. Very recently, there has been some the- 
oretical work concerning Feshbach resonances in BECs 
in optical lattices [211 122 . [23| . This situation is also the 
subject of current experimental investigations 

In the present paper, we use an averaged CP equa- 
tion [2^ I23 to examine gap solitons in BECs trapped in 
optical lattices and under the influence of Feshbach reso- 
nance management. Using an envelope- wave approxima- 
tion, we derive coupled-mode equations describing the 
slow dynamics of gap solitons. We provide an analyt- 
ical construction of gap soliton solutions and examine 
their spectral stability with numerical computations of 
eigenvalues. We then describe the time-evolution of gap 
solitons with the averaged and original CP equations. 
Finally, we summarize our results. 



THE AVERAGED GROSS-PITAEVSKY 
EQUATION 

We consider the non-dimensional CP equation for 
trapped BECs under Feshbach resonance management. 



and the tildes have been dropped from Eq. j^J. The 
nonlinear coefficient can be written 1201 



where the normalized independent variables are 
Jimx ~ t 



(2) 
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and the potential for optical lattices is [27 



V{x) = 2eVo cos{ljx) 



(3) 



(4) 



Here, 70 is the mean value of the nonlinearity coefficient; 
7(t), with r = t/e, is a mean-zero periodic function with 
unit period; e <C 1 is a small parameter describing the 
strength of the Feshbach resonance management; w is 
the wavenumber of the optical lattice; and eVo is a small 
parameter describing the strength (amplitude) of the lat- 
tice. Because ^ changes rapidly in time (i.e., the man- 
agement is strong), it is reasonable use a model without 
dissipation [31 . We note that the two small parameters e 
and e in Eq. Q are due to two different physical sources 
and can be uncorrelated. 

One can exert very precise control over optical lat- 
tice strengths and wavenumbers experimentally [l^. In 
particular, both strong and weak lattices can be imple- 
mented; we will consider a weak optical lattice, so that 
e in Q is small. In experiments, the scattering length a 
can also be adjusted either nonadiabatically (strong non- 
linearity management), which is the situation discussed 
in the present work, or adiabatically (weak nonlinear- 
ity management). In general, different dynamics can 
occur dependin g o n how rapidly a is adjusted, as dis- 
cussed, e.g., in j29|. There have already been some ex- 
periments in which Feshbach resonances are applied to 
BECs in the presence of optical lattices jl^l and others 
are being planned by multiple experimental groups. As 
mentioned above, a conservative model is appropriate for 
strong nonlinearity management. For weak nonlinearity 
management, however, the CP equation ^ needs to be 
augmented by dissipation terms, as three-body recombi- 
nation leads to experimental losses [HI E3|. Therefore, 
we will only consider strong management, so that e in 
© is small. 

The small parameter e can be used to simplify the 
time-periodic CP equation ^ with an averaging method. 
Using the time-averaging procedure from [2!tI |. one can 
look for an asymptotic solution to the CP equation (0 
of the form 



V'(a;,t) =e*"^-i(^)l"''(=^'*) [m(x, t) + 0(e)] 



(5) 



where u is the complex- valued amplitude and 7_i(t) is 
the mean-zero anti-derivative of "fir). Within a regular 
averaging procedure (see the review in [3), we obtain 
the averaged CP equation, 

iut = —Uxx + 2eVb cos{lljx)u + 7o|upu 

-l![{{\u\')xr + 2\u\\\uWx]u, (6) 
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where 71 is the standard deviation of the nonhnearity 
coefficient. The model ^ provides a starting point for 
our analysis of gap solitons in BECs in optical lattices 
under Feshbach resonance management. Without such 
management (71 = 0}, gap solitons of the GP equation 
® were studied in [23 using Floquet theory, multiple- 
scale expansions, beyond-all-orders theory, and Evans- 
function computations of eigenvalues. We will mainly 
consider the opposite case when 71 7^ but 70 = 0. 



COUPLED-MODE EQUATIONS 

We are interested in modeling gap solitons supported 
by subharmonic resonances between the periodic poten- 
tial and spatio-temporal solutions of the averaged GP 
equation ©. To simplify the model, we use the second 
small parameter e and obtain coupled- mode equations, 
which averages the spatially periodic nonlinear equation 
near a spectral gap of its associated linearization, 



nian, with symmetric potential energy 



iut = —Uxx + 2eVo cob{ljjx)u . 



(7) 



In the limit of small e, the spectral gaps all become nar- 
row and the first spectral gap occurs at first order in 
e. Using the space-averaging technique of Ref. [s^l, one 
can look for an asymptotic solution of the averaged GP 
equation © in the two-wave form 

u{x,t) = (yl(X,r)e*'^°^-''^«* 

+B(X,r)e-*'^o^"''^«* + 0(e)) , (8) 

where A and B are complex-valued amplitudes, X = 
ex and T — et are slow variables, and to = 2clIo- This 
wavenumber ratio indicates that we are studying 2 : 1 
subharmonic resonances [e.g. the first spectral gap of 
equation Using the regular asymptotic procedure 

from fs^l , we obtain a system of coupled-mode equations 

i{At + ujAx) = VoB + 7o(|A|2 ^ 2|B|2)A 

+ 27iV(2|^|' + |Sn|BpA, 

i{BT - LoBx) = VoA + 7o(2|A|2 ^ j^p^^ 

-F27iV(|A|2+2|B|2)|A|2b. (9) 

The first equation governs the left-propagating wave 
A, whereas the second equation governs the right- 
propagating wave B. The two waves interact with a cubic 
cross-phase modulation from the mean value of the scat- 
tering length and with a quintic cross-phase modulation 
from the standard deviation of the scattering length. The 
latter effect represents the main contribution of Feshbach 
resonance management in trapped BECs. 

The system of coupled-mode equations © is Hamilto- 



W{A,B) 



Vo{AB 
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AB) 
{\A\^ + 4\Af\B\ 



\B[ 



+ 2YiUj'\A\'\B\'{\A\' + \B\'). 



(10) 



Additionally, Eq. (|10|l satisfies the assumption on sym- 
metric potential functions used recently for analyzing the 
existence and stability of gap solitons While the pre- 
vious work concerned gap solitons in coupled- mode equa- 
tions with cubic nonlinearity, we will focus on the new 
effects that arise from the quintic nonlinear terms. These 
effects correspond to mean-zero Feshbach resonance man- 
agement, which affects the propagation of gap solitons in 
optical lattices. We can see that the last term of W 
in is positive definite, similar to the second term 
with 7o > 0. Therefore, Feshbach resonance manage- 
ment leads to defocusing effects on the propagation of 
gap solitons in periodic potentials. The defocusing role of 
Feshbach resonance management was studied recently in 
[32| | in the context of blow-up arrest in multi-dimensional 
GP equations. 

One determines the linear spectrum of the coupled- 
mode equations ijHl from the linearized system in Fourier 
form, {A,B) - e^-^-^^^"^, where ft = ±yJV^ ^ iJ^K"^ . 
The spectral gap exists for |ri| < |Vo| and corresponds to 
the first spectral gap associated with the periodic poten- 
tial in Eq. Q. The lower (upper) spectral band of the 
coupled- mode equations Q for < — |Vb| (for > |Vo|) 
corresponds to the first (second) spectral band of the pe- 
riodic potential in (0. 

The coupled-mode system © can be reduced to a non- 
linear Schrodinger (NLS) equation near the band edges 
of the linear spectrum. Using the asymptotic represen- 
tation 



ZVq 



(11) 
, (12) 



where S, = fiX, ^ = /i^T, and fi is a small parameter 
for the distance of ft from the band edges ±|Vo|, one can 
reduce the coupled-mode equations © with 70 = to 
the quintic NLS equation: 



6-ffLu^\W\^W . 



(13) 



The quintic NLS equation p3(l is focusing near the lower 
spectral band with = — | Vq | and is defocusing near the 
upper spectral band with ft = \Vo\. Therefore, the gap 
soliton solutions bifurcate from the lower spectral band 
via a local (small-amplitude) bifurcation. They termi- 
nate at the upper spectral band via a nonlocal (large- 
amplitude) bifurcation, similar to gap solitons in the GP 
equation with a periodic potential |27l] . 



4 



The derivation of the quintic NLS e quat ion H13() con- 
firms the predictions of the recent paper |33| on the power 
threshold for one-dimensional gap solitons in the case 
7o = and 71 ^ 0. The existence of the power threshold 
near the local bifurcation limit was computed numerically 
in Ref . js^ . Because the quintic NLS equation exhibits a 
similar power threshold for NLS solitons (see, e.g., [l^l), 
the numerical fact is now confirmed from the perspective 
of asymptotic theory. We note that the coupled-mode 
equations (jnj with 70 7^ reduce to the cubic NLS equa- 
tion, which does not exhibit the power threshold for NLS 
solitons. 



GAP SOLITONS 

We simplify the construction of exact gap soliton solu- 
tions to the coupled-mode equations (j^ by normalizing 
Vo = — l,iu; = l (a standard scaling transformation can 
be employed for this purpose) and defining = 27^0;^. 
We construct gap soliton solutions by separating vari- 
ables into time-periodic and spatially localized solutions 
of the coupled- mode equations ((HJ: 



A{X,T) = a(X)e 



B{X,T) ^ b{X)e 



iUT 



Because of the symmetry in the potential function, 
W{A,B) = W(B,A). the gap soliton solutions satisfy 
the constraint b = a [SlJ, so that a{X) solves the follow- 
ing nonlinear ordinary differential equation: 

ia' + fla + d — 37o|a|^a -f 3(T'^|a|*a . 

Converting the function a{X) to polar coordinates, 



(14) 



a{X) =. v/QW exp[-ie(X)/2] , 
we obtain the second-order system, 
Q' = --2gsine, 

e' = -2r2-2cose-|-67oQ + 6(T2Q2 (15) 
This system has the first integral 



-QQ - cos e + -joQ^ +a^Q'^., 



(16) 



where E = from the zero boundary conditions Q{X) — > 
as \X\ — + 00. 

In the remainder of this paper, we consider the case 
7o = and 71 7^ 0, which shows the effects of Feshbach 
resonance management on the existence and stability of 
gap solitons. For the case 70 7^ and 71 = 0, exact 
analytical solutions for gap solitons are available j35j and 
the stability problem has been analyzed numerically (see 
[si'l). When 70 — 0, the second-order system lfT3)l reduces 
to the first-order differential equation. 



e' = 4(f7-|-cose) 



(17) 



The function Q{X) is found from QiX) with the relation 



i} + cos Q 

Q = — — 



> 0. 



(18) 



Using the technique from Appendix A of [31| , we find the 
exact analytical solution of Eq. H17() and obtain 

cosh^ (2/3x)~7sinh' {2Px) 
cos 6 = . T , . . T , , (19) 



where 



cosh^ (2/?x) -t-7sinh2 (2/3a;) ' 



and |n| < 1. Substituting fH^ into ifT^ then gives 

1 1 + n 



Q - 

o"^ cosh^ {2i3x) -I- 7 sinh^ (2/3a;) 



so that 



a{X) 



^7(cosh (4/3 a:) -17) 
V^[cosh(2/3X) + i/7sinh(2/3A:)] 



and 



\a? = 



/3 



cr^cosh(4/3A:) - n 



(20) 



(21) 



(22) 



The limit f2 — > — 1 yields the small-amplitude gap soliton 
|ap sech(2/3A'), which satisfies the focusing 

quintic NLS equation (|13|l . At this local bifurcation limit, 
the power of the gap soliton has a threshold. 



P 



\Bn dx ^ Po 



such that the power P is bounded from below by 
the limiting value Pq. The opposite limit Vl +1 
yields the large-amplitude (singular) gap soliton |ap — > 
(/3/cr-\/2)csch(2/3Ar), which satisfies the corresponding de- 
focusing quintic NLS equation (|13|) . Thus, in accor- 
dance with the asymptotic reduction to the quintic NLS 
equation (|13|) . the family of gap soliton solutions of the 
coupled-mode system (j^l bifurcates from the lower spec- 
tral band (fJ = —1) and terminates at the upper spectral 
band [VL = +1). 



Stability 

The spectral stability of the gap soliton 1)21(1 follows 
from the linearization 



iO,T 



A{X,T) 

A{X,T) 
B{X, T) = e 

B{X,T) = e'''^'(a{X) 



a{X) 
e— [a{X)- 
{d{X) 



hC/i(X)e^^) 
C/2(X)e^^) 
HC/3(^)e^^) 
- C/4(X)e/^) 



(23) 



5 



The vector U = (C/i, t/2, C^3, Kt)"^ solves the hnear eigen- 
value problem, 



H^XJ = iXsV, 



(24) 



where s = diag(l, — 1, 1, — 1) is a diagonal matrix. The 
linearized energy operator has the form 



D{dx) + ViX). 



where 



D = 



and 



/ -n- idx 


-1 





-1 







-1 





-SI + idx 

-n + idx 
-1 -n-idx J 



(25) 



/ 5|o|4 2|a|2a2 A\a\''a' 4|a|4 \ 

2|a|2a2 5|a|4 4|a|4 A\a\''a^ 

4|a|2a2 4|o|4 5|a|4 2|a|2a2 

V 4|a|4 4|a|2a2 2|a|2a2 5|a|4 / 



Using the block-diagonalization method from [3l|, we 
employ the orthogonal similarity matrix 



\/2 



/ 1 1 \ 
10 1 
10-1 

V 1 -1 / 



that simultaneously block-diagonalizes the energy oper- 
ator Hi^, 



S HuiS — 



H+ 
iJ_ 



and the linearized operator sH^ 



iJ_ 



H+ 

where -ff± are two-by-two Dirac operators 



iL . 



(26) 



(27) 



i/4 



-n-idx + 9(r^\a\ 
6a^\a\^a^-l 



—fl — idx 

1 -2a2|a|2a2 



6fT2|a|2a2_i 



1 - 2cr2|a|2a2 



(28) 



(29) 



Eigenvalues of the operators L, H+, and H- are detected 
numerically with the Chebyshev interpolation method 
[sH - The main advantage of the Chebyshev grid is that 
clustering of the grid points occurs near the end points 
of the interval. This prevents the appearance of spurious 
complex eigenvalues that may otherwise arise from the 
discretization of the continuous spectrum. Moreover, by 
using the block-diagonalization in H26|) - (|27|l . we are able 



to reduce the memory constraints and double the speed 
of the numerical computations (see the details in |3l|). 

The numerical eigenvalues of the operators L, 
and are displayed in Fig. ^ for six different values 
of the parameter fi. When fl is close to the local bi- 
furcation threshold (e.g., for ft = —0.9), the operator 
L has a four-dimensional kernel at A = and a pair of 
small purely imaginary eigenvalues near A = 0. The pair 
of purely imaginary eigenvalues originates from the six- 
dimensional kernel of the linearized quintic NLS equation 
((T^ (see, e.g., Ref. 0). In this case, the operator if+ 
has no isolated non-zero eigenvalues, whereas the oper- 
ator H- has a simple isolated non-zero eigenvalue. The 
eigenvalue of still exists at = —0.7, but it disap- 
pears before — —0.3 because it collides with the end 
point of the continuous spectrum of The pair of 

eigenvalues of L survives at fl = —0.3, but it disappears 
at SI = because it collides with the end points of the 
continuous spectrum of L. When S7 > 0, the pair of com- 
plex eigenvalues bifurcates in the spectrum of L from the 
end points of the continuous spectrum of L. The complex 
eigenvalues of L bifurcate simultaneously with a simple 
isolated non-zero eigenvalue of the operator H^. The 
pair of complex eigenvalues of L and the isolated non- 
zero eigenvalue of H- persist for larger values of SI (e.g., 
for S7 = 0.3 and SI = 0.7). The eigenvalues just discussed 
are labeled in the figure as I and II. 

In sum, the gap solitons of the coupled-mode system 
Q with 7o = are spectrally stable for < and spec- 
trally unstable due to complex unstable eigenvalues for 
S7 > 0. This behavior is similar to the stability analysis 
of the coupled-mode system Q with 70 7^ and 71 = 
(see [33 and references therein). 



NUMERICAL SIMULATIONS 

To confirm the accuracy of the coupled-mode the- 
ory, we corroborate the instability of gap solitons for 
S7 > predicted from the analysis of the system ^ with 
full numerical simulations of the averaged GP equation 
© . We consider the trigonometric management function 
7(t) ~ cos(27rT) and fix the other parameters as follows: 
70 - 0, 71 = I/V2, Vo = -1, cc; = 1, and e = 0.1. 
We vary the parameter S7 for the gap soliton solutions, 
focusing on computations with S7 — —0.5 and S7 — 0.5. 
According to analysis of the coupled- mode equations 
the value = —0.5 corresponds to the case of stable gap 
solitons, and the value S7 = 0.5 corresponds to the case 
of unstable gap solitons. The initial condition for all sim- 
ulations is selected to be a perturbation of the leading- 
order two-wave approximation |(SJ) with the gap soliton 
solutions 

We integrated the averaged GP equation © using a 
finite-difference approximation (with 480 grid points) in 
space and a Runge-Kutta integration scheme (with step 
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size h = 0.001) in time. The perturbations of the initial 
gap soUtons are of the same functional form as the gap 
solitons but with larger amplitudes. 

Figure El shows the time-evolution of a stable gap soli- 
ton with f2 = —0.5. The stationary gap soliton persists 
in the full dynamics of the averaged GP equation 10, 
in agreement with the stability analysis of the coupled- 
mode system 

Figure |2| illustrates the dynamics of an unstable gap 
soliton with O = 0.5. We observe an asymmetric beat- 
ing between different localized waveforms. The localized 
wave is not destroyed in the unstable case, but rather 
undergoes shape distortions due to the oscillatory insta- 
bility. This behavior agrees with the stability analysis 
in the coupled-mode system (jSl, as the unstable eigen- 
values of gap solitons with 17 > are complex- valued 
(see Fig. Pl. While the perturbation used in Fig. |31 is 
large, note that a perturbation with the same percentage 
difference in the soliton amplitude does not lead to an 
instability in Fig. [3 

We also examined the dynamics of the solutions con- 
structed using coupled-mode theory under full numerical 
simulations of the original GP equation (|2J). In these 
simulations, which were also performed using a finite- 
difference approximation (with 480 grid points) in space 
and a Runge-Kutta integration scheme (with step size 
h — 0.001) in time, the external potential includes con- 
tributions from a small harmonic trap in addition to the 
optical lattice in Q. Hence, the potential in these sim- 
ulations is given by V{x) — V{x) + VhX^ ■ The initial 
conditions and parameters are the same as for the simu- 
lations of the averaged GP equation, with the additional 
values £ = 0.25 and Vh = 0-01. In Fig. ^ we show the 
spatio-temporal evolution of the constructed gap solitons 
under equation (0). As expected, the solitons that we 
found to be stable persist longer with respect to the full 
GP dynamics than those determined to be unstable. 



furcation threshold). We then showed with numerical 
simulations of the averaged Gross-Pitaevsky (GP) equa- 
tion that unstable gap solitons exhibit beating between 
different localized shapes, thereby confirming the stabil- 
ity results predicted from the coupled-mode theory. We 
corroborated this further with numerical simulations of 
the original GP equation, which show that the stable gap 
solitons persist much longer than the unstable ones. 

We note that gap solitons in the GP equation with a 
periodic potential © and the coupled-mode system © in 
the case of no Feshbach resonance management have been 
studied recently in [i^l and , respectively. We can see 
from comparing the previous and new results that Fes- 
hbach resonance management leads to a new effect with 
respect to the existence of the power threshold near the 
local bifurcation limit ^3 • On the other hand, there are 
not many differences in the stability results, as Feshbach 
resonance management does not stabilize gap solitons far 
from the threshold of local bifurcations (which are known 
to be unstable 0|). We have also confirmed that Fes- 
hbach resonance management leads to defocusing effects 
on the propagation of gap solitons, which is relevant for 
blow-up arrest in multi-dimensional problems Is^ . 
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CONCLUSIONS 



In conclusion, we studied Feshbach resonance man- 
agement for gap solitons in Bose-Einstein condensates 
trapped in optical lattice potentials. We applied an 
envelope wave approximation to the averaged Gross- 
Pitaevsky equation with a periodic potential to yield 
coupled-mode equations describing the slow BEG dynam- 
ics in the first spectral gap of the optical lattice. We de- 
rived exact analytical expressions for gap solitons with 
zero mean scattering length. In this situation, Feshbach 
resonances are employed to tune a condensate between 
the repulsive and attractive regimes (corresponding to 
their usual experimental application). Applying Gheby- 
shev interpolation to the coupled-mode equations, we 
showed that these gap solitons are unstable above the 
center of the first spectral gap (far from the local bi- 
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FIG. 1: (Color online) Eigenvalues and instability bifurcations 
for the operators L, H+, and H- in ^ and The 

parameter values are Q = —0.9 (top left), Q — —0.7 (top 
right), n = -0.3 (middle left), Q = 0.0 (middle right), Q = 
0.3 (bottom left), and Q. = 0.7 (bottom right). 




FIG. 2: (Color online) Stable time-evolution of a gap soli- 
ton with = —0.5 in the averaged GP equation (Left) 
Spatio-temporal evolution of \u{x,t)f. (Right) Spatial pro- 
files of \u{x,t)\^ for t = 0, t = 200, t = 400, and t = 800. 




FIG. 3: (Color online) Unstable time-evolution of a gap soli- 
ton with n = 0.5 in the averaged GP equation (Left) 
Spatio-temporal evolution of |u(a;,i)p. (Right) Spatial pro- 
files of \u{x,t)f for t = 0, t = 200, t = 400, and t = 800. 
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FIG. 4: (Color online) Spatio-temporal evolution of \ip{x,t)f , 
with initial conditions given by the constructed gap solitons, 
under the dynamics of the full GP equation ||5J. (Left) Q = 
-0.5. (Right) n = 0.5. 



